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Abstract 

The field mixing that manifests broken particle-hole symmetry, is studied for a 2-D 
asymmetric lattice gas model having tunable field mixing properties. Monte Carlo simu- 
lations within the grand canonical ensemble are used to obtain the critical density distri- 
bution for different degrees of particle-hole asymmetry. Except in the special case when 
this asymmetry vanishes, the density distributions exhibit an antisymmetric correction 
to the limiting scale-invariant form. The presence of this correction reflects the mixing 
of the critical energy density into the ordering operator. Its functional form is found to 
be in excellent agreement with that predicted by the mixed-field finite-size-scaling theory 
of Bruce and Wilding. A computational procedure for measuring the significant field 
mixing parameter is also described, and its accuracy gauged by comparing the results 
with exact values obtained analytically. 

PACS numbers 64.70, 64.70F 



1 Introduction 



It has long been appreciated |IJ that the lack of symmetry between the liquid and vapour 
phases of a fluid leads to a mixing of the temperature and chemical potential in the two 
relevant scaling fields close to the critical point. This is in marked contrast to the situation 
for models of the Ising symmetry such as the ordinary lattice gas ||, which as a consequence 
of their so-called 'particle-hole' symmetry^ exhibit no field mixing. Although the reduced 
symmetry of fluids has no consequences for their universal properties (which for systems with 
short-ranged interactions correspond to the Ising universality class), it is predicted that certain 
non-universal effects are engendered by field mixing. Principal among these, is the existence of 
a weak energy-like singularity of the coexistence diameter on the approach to criticality. The 
presence of this singularity, now firmly established experimentally || |4[] , constitutes a failure 
for the century-old empirical 'law of rectilinear diameter'. Indeed, some theoretical progress 
has been made towards an understanding of the microscopic factors governing the amplitude 
of the diameter singularity || || . 

Recently however, a new finite-size-scaling theory has been developed that relates the mixed 
character of the fluid scaling fields to the interplay of the near-critical energy and density 
fluctuations |7|, [|. The theory predicts that as a result of field mixing, the energy operator 
features in the critical density distribution, giving rise (as outlined below) to a correction to the 
limiting (large L) universal form of the density distribution. This correction is subdominant 
to the limiting form and dies away with increasing L. To leading order in the theory, its 
functional form is prescribed by independently known functions characteristic of the Ising 
universality class. Moreover the symmetry of the correction (an antisymmetric function) differs 
from that of the limiting form (a symmetric function). Its presence in the density distribution 
is therefore a potentially distinctive signature of the field mixing phenomenon, one that can 
in principle be isolated and analysed by means of computer simulation measurements of the 
density fluctuations. Indeed, the potential utility of the theory in facilitating simulation studies 
of field mixing was clearly demonstrated in an extensive Monte Carlo investigation of the 2-D 
Lennard- Jones fluid near the liquid- vapour critical point ||. Measurements of the near-critical 
density distribution yielded an antisymmetric field mixing component that mapped quite well 
onto the predicted universal form. In addition, a computational prescription was set out for 
estimating the more significant of the two field mixing parameters (that which controls the 
extent to which the chemical potential features in the thermal scaling field). 

Notwithstanding the successes of the study reported in |7|, [|, one is still confronted with 
a number of difficulties if one wishes to assess the quantitative validity of the new theory. 
The most severe problem is the computational burden imposed by simulations of realistic fluid 
models (such as the Lennard- Jones system), which owing to their continuous, long-ranged 
interaction potential, require intensive floating-point calculations. The computational demands 
of such simulations far exceed those of lattice-based particle or spin systems. Moreover, since 
exact values of the critical couplings are not generally available for realistic fluid models, 
it is difficult to probe the asymptotic critical region. This can complicate the isolation of 
the field mixing correction to the density distribution, which must be identified amidst other 
corrections associated with small departures from criticality. The analytical intractability 
of realistic fluid models also precludes an assessment of the accuracy of any computational 

*A lattice gas possesses particle-hole symmetry if the Hamiltonian, a function of the site occupation numbers 
tJi = 0, 1 satisfies the relation 7i — > H+A^i <Ji+B (with A, B constants) under the transformation g l — > 1 — er^. 
Real (continuous) fluids lack particle-hole symmetry ipso-facto since a hole is not defined. 
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procedure for measuring the significant mixing parameter. 

Clearly, in order to facilitate a more detailed assessment of the mixed-field finite-size-scaling 
theory, it would be beneficial to work with a model system having fewer of the drawbacks 
listed above. One such system is the 2-D asymmetric lattice gas model, originally proposed 
by Mermin, which whilst being comparatively much easier to tackle computationally than 
e.g the Lennard- Jones fluid, is also analytically solvable for the critical point couplings and 
field mixing parameters. It therefore provides an ideal test-bed for the theory, permitting 
an accurate determination of the field mixing correction, and providing a benchmark against 
which, the accuracy of any computational method for determining the field mixing parameters 
can be gauged. 

Below we describe the results of field mixing studies of Mermin's asymmetric lattice gas 
model. The layout of the paper is as follows: the model is described in section |2| and set within 
the framework of the mixed-field finite-size-scaling theory. Exact values for the significant 
mixing parameter of the model are also calculated. In section |^, Monte Carlo measurements 
of the critical density distribution are presented. The antisymmetric component of these 
distributions is isolated and compared with the predicted form of the field mixing correction. 
Estimates of the significant mixing parameter are also deduced from the simulation data and 
compared with the exact results. Section f| details our conclusions. 

2 Background 

The asymmetric lattice gas model forming the focus of the present work, was first proposed by 
Mermin as an example of a system exhibiting a singular coexistence diameter^. The model 
consists of an ordinary 2-D square lattice gas (whose nearest neighbour coupling we denote J) 
in which particle-hole symmetry is manually destroyed by forbidding occupation of sites whose 
row and column numbers are both even. By this action one creates two types of sites: odd- 
odd sites having coordination number 4, and odd-even (or even-odd) sites having coordination 
number 2. It is straightforward to show that the particle- hole symmetry that obtains in the 
ordinary lattice gas is equivalent to the requirement that all sites have the same average energy 
environment. The presence of two sets of inequivalent sublattices in the asymmetric lattice 
gas, clearly violates this condition and leads to field mixing. It transpires, however, that if one 
introduces an additional coupling K, between atoms on the sublattice of odd-odd sites, then 
the degree of particle- hole asymmetry can be tuned. Indeed for the special choice K = — J/2, 
the average energy per site becomes equal for both sublattices and consequently particle-hole 
symmetry is once more restored. 

Aside from its field mixing properties, the chief asset of Mermin's asymmetric lattice gas, is 
its analytic tractability. The grand partition function of the asymmetric model can be related 
by means of analytic transforms to that of the ordinary lattice gas, for which in turn a wealth 
of exact results are known in two dimensions 0. Specifically, one finds ||: 

£l(lM,T) = (l + e^^UfcT) (1) 

where Q is the partition function of the asymmetric model and \x and T are the chemical 
potential and temperature respectively. Bars denotes quantities in the ordinary lattice gas and 

2 We note in passing, that the model also possesses a rich phase structure that has been investigated in 
detail by other workers £u], 0, 0, H| 
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N = N/3 where N is the number of allowed sites in the asymmetric model. 

Introducing the dimensionless chemical potential £ = fx/ kbT and dimensionless coupling 
parameters 77 = J / k b T, A = Kj fc b T, equation IT] leads to the following relationships M : 



£ 
V 



£ + 41n 
A + ln 



1 + e 



l + e« 
l + eg)(l + eg +2 ") 
(1 + e?+") 2 



(2a) 
(2b) 



where £ = Ji/kbT and 77 = J /kbT are respectively the dimensionless chemical potential and 
dimensionless nearest neighbour coupling constant of the ordinary lattice gas. 

Now, it transpires that in the ordinary lattice gas the liquid-vapour coexistence line is 
specified by £ = — 2rj, while the critical point that terminates this line is given by the solutions 
to the relations: 



Setting £ = - 
metric model: 



sinh(r7 c /2) = 1, £ c = -2rJ c (3) 
■27} in equations ^a| and |2b] then yields the coexistence condition for the asym- 



£ + 2A + 2 In 



1 + e 







(4) 



which represents a surface in the space of £, 77 and A. Note however, that since A and 77 both 
enter only as multiplicative factors in the configurational energy (see equation^ below), this 
coexistence surface can be represented in terms of a family of coexistence curves in the spacef] 
of £ and rj, each curve being parameterised by a different value of the coupling ratio A/77. 
This ratio is consequently the crucial parameter controlling the coexistence and field mixing 
properties. 

We shall also find it useful to obtain the critical density of the asymmetric lattice gas, 
which can be calculated by appeal to the relation: 



Pc 



1 dlntt(n,T) 
N ' 



from whence, one obtains 



Pc 



l + e 




(5) 



(6) 



where the subscript c on the derivatives signifies that they are to be evaluated at criticality. 
The quantities u c and p c are respectively the critical energy density and number density of the 
ordinary lattice gas model, the particle-hole symmetry of which implies ~p c = 0.5. The value of 
u c is also known exactly; it is related to the critical energy density of the 2-D Ising model 
by the relation u c = |[2 + where from Onsager's solution, u\ = y/2/2. The derivatives in 
equation Fjl can be calculated straightforwardly to yield 



3 Note that we might equally well have chosen A instead of rj as the independent coupling variable. 
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d£ 1 + 5et +r > - 3e e + e 2 ^ 

dl ~ (l + e5)(l + e^) ( } 

drj e i _ e 2?+r? + 2e 2 ^ + e t +2r > - e 2 ^ - 2e«+" 

d£ ~ (1 + e$)(l + et+ 2 i)(l + et+i) ^ ' 

Substituting for TJ C and £ c in equations |2a] and and feeding the results for r\ c and £ c into 
equation || one readily obtains the critical density p c as a function of A/77. This relationship is 
shown in figure [I] for values of A/77 in the range (—1/2, 1/2), which, as will be seen, encompasses 
a wide range of field mixing behaviour. The figure clearly demonstrates that when A/77 = —1/2 
the critical density is that of the ordinary lattice gas (p c = 0.5), thus confirming that particle 
hole symmetry is restored for this value of the coupling ratio. It is also evident that increasing 
A/77 causes the critical density to decrease monotonically. The nature of this density shift finds 
illustration in the simulation results to be described later. 

We turn now to the field mixing properties of the asymmetric lattice gas model, which we 
analyse within the framework of the mixed-field finite-size-scaling theory of references 0, ||. To 
this end, we consider a 2-D system of side L, having a maximum available volume V = 3L 2 /4. 
The configurational energy $({<r}) is given by 

= Wi^j + Xa niVn (8) 

<hj> [m,n] 

with <ji = 0,1. The site indices % and j are taken to run over all allowed nearest neighbour 
sites, while m and n run only over nearest neighbours on the sublattice of odd-odd sites. We 
assume also that the system is thermodynamically open so that the particle number density 
p = J2i a i/V can fluctuate. In this paper, we shall be concerned with the statistical behaviour 
of both the number density p, and the configurational energy density u = 1 /_1 ?7 _1( I ) ({(t}), the 
latter of which we write in units of the coupling parameter 77. 

For a given choice of the coupling ratio A/77, the critical point is located by critical values 
of the reduced chemical potential £ c and reduced coupling rj c . Deviations of £ and r\ from 



their critical values determine the size of the two relevant scaling fields [H|. In the absence of 
particle-hole symmetry, it is expected [0] that the relevant scaling fields comprise (asymptoti- 
cally) linear combinations of these deviations: 

r = Vc ~ V + s(f - 6c) h = £-£ c + r(r] c -ri) (9) 

where r is the temperature-like scaling field and h is the field-like scaling field. The parameters 
s and r are system-specific quantities controlling the degree of field mixing. In particular, 
r is identifiable as the limiting critical slope of the coexistence curve in the space of £ and 77 
[0]. The role of the parameter s (which controls the degree to which the chemical potential 
features in the thermal scaling field) is, however, more significant: it determines the size of the 
diameter singularity. We term s the significant mixing parameter. 

Conjugate to the two relevant scaling fields are scaling operators £ and A4, which are found 
to comprise linear combinations of the energy and number densities ||: 

£ = j^gf [u — rp] M = j^gf [p — su] (10) 
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where £ is the energy-like operator and Ai the ordering operator. In the Ising context (for 
which s = r = 0), Ai is simply the magnetisation while £ is the energy density. 

Near criticality, the joint probability distribution Pm,£(M-,£) of the operators £ and Ai is 
expected to exhibit scaling behaviour. In particular, in the limit of large system size L, the 
distribution of the ordering operator Pm{M) = J PM,s(M,£)d£ is expected to be describable 
by a finite-size-scaling relation having the form [ |15|j : 

p M (M) ~ a M - x L d - XM p M {a M - 1 L d -^bM, a M L XM h, a £ L X£ r) (11) 

where 5Ai = Ai — Ai c and the function Pm is predicted to be universal, modulo the choice for 
the scale-factors and a £ of the two relevant fields, whose scaling indices are Am = d — (3/v 
and Ag = 1/v respectively. This scaling form (which has its basis in the renormalisation 
group scaling properties of the multi-point correlation functions [[U|) is well-supported by 
Monte Carlo studies of 2-D Ising and 4> 4 models, where the ordering operator Ai is simply the 
magnetisation p"7 . 



In asymmetric systems, the mixing of the energy density into Ai precludes direct simulation 
measurements oipM{Ai), since in general the value of s will not be known a-priori. Instead, 
it is expedient to focus on the distribution of the density, which is obtained from the joint 
distribution of the mixed operators by integrating over the energy spectrum and expanding in 
the mixing parameter s : 

Pl(p) = J pmAp - s£ i £ ) d£ ( 12 ) 

from which one finds M 

p L (p) ~ a A1 - 1 L d - A ^(a A4 - 1 L d - A ^[p - p c ],a M L x ^h,a £ L X£ T) + Ap L (p) (13) 

with 

d 

^Pl(p) = -8 j- {p L (p) [< u{p) > -u c - r{p - p c ))} + 0(s 2 ) (14) 

where < u(p) > (hereafter referred to as the energy function) is the mean energy density for 
a given p. 

The function Ap^(p) describes (to linear order in s) the component of the critical density 
distribution associated with field mixing. Precisely at criticality, equation [13| may be written 
in the form: 



Pl(p) ^ a M - x L^ 



PU(*) - sa £ a M - x L-^-^^- (p* M (x)?{x)) + 0{s 2 ) 

(15) 

where p* M (x) = Pm( x i 0, 0) and 



d\np M (x,0,z) 
e (x)= — 



dz 



= a e - x L d -^\< u{p) > -u c - r{p - p c )\ + 0{s) (16) 

2=0 

is a universal function whose form has (like that of p^) been previously established in Monte 
Carlo studies of the critical 2-D Ising model where it is simply the energy function for 
the magnetisation. 
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To leading order in s, the critical density distribution can thus be expressed as a sum of two 
independently-known universal components. The first of these, p* M (x), is a function having 
the same form as the critical magnetisation distribution of the Ising model. The second, 

(p* M (x)e*(x)), is a function characterising the critical energy operator and represents (to 
linear order in s) the field mixing contribution to the density distribution. This field mixing 
term is down on the first term by a factor T J -( 1 - a -i 3 )/ l/ anc l therefore represents a correction 
to the large L limiting behaviour. Given further the symmetries of e(x) and p* M (x), both of 
which are even (symmetric) in the scaling variable x, the field mixing correction is the leading 
antisymmetric contribution to the density distribution. Accordingly, it can be isolated from 
measurements of the critical density distribution, simply by antisymmeterising around p c . 



In addition to furnishing the functional form of the field mixing correction, equation [14 
also constitutes a computational prescription for estimating the significant mixing parameter 
s . The value of s is simply the single scale factor required to match the critical point form 
of the measured function Ap L (p) = —s-^ {pl(p) [< u(p) > —u c — r(p — p c )}} to the measured 
antisymmetric component of the critical density distribution. In the present case, the accuracy 
of this procedure can be gauged, since the exact value of s is obtainable analytically. 

An exact calculation of s follows from the observation that in the ordinary lattice gas, the 
field-like scaling field h coincides with the line fj = fj c in the space of £ and TJ. It follows that 
in the asymmetric model, the direction of h in the space of £ and rj can be obtained from 
equation |2| by setting fj = fj c , A = \(r)) and solving for r\. The value of s is then given by 

(17) 

where the derivative is to be evaluated at criticality. The results of this calculation are displayed 
in figure |2] as a function of A/77 in the range (—1/2,1/2). As anticipated, s vanishes for 
A/77 = —1/2. We note also the presence of a broad minimum in the value of s at A/77 — — 0.08. 




3 Monte Carlo studies 
3.1 Computational details 

The Monte Carlo simulations reported here, were all performed using a Metropolis algorithm 
within the grand canonical ensemble. Two system sizes were studied, having linear dimension 
L = 20 and L = 30. Periodic boundary conditions were employed throughout. The basic 
observables recorded were the probability distribution of the number density Pl{p) and the 
energy function < u(p) >. The distribution of the density was obtained initially as a histogram. 
The energy function was accumulated as an average of the energy for each value of p explored 
in the course of the simulation. All simulations were performed at the exact critical point of 
the model, obtained (for a given choice of A/77) front equations |2a|, ^b] and ^. The L = 20 
systems consisted of 1 x 10 6 lattice sweeps for equilibrium, followed by a sequence of 2 x 10 5 
observations with 100 sweeps between each observation. For the L = 30 system, equilibration 
times of 4 x 10 6 sweeps were used with 2 x 10 5 observations separated by 200 lattice sweeps. In 
each instance the whole procedure was repeated 12 times to test the statistical independence 
of the data and to assign statistical errors to the results. 
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3.2 Results 



Measurements of the critical density distribution were collected for the L = 20 system at five 
distinct values of A/77, namely —1/2, —1/4,0, 1/4, 1/2. For the L = 30 system, measurements 
were made with A/77 = —5/12 and A/77 = 0. The L = 20 distributions are shown in figure [3], 
where each has been normalised to unit integrated weight. With the exception of the case A/77 = 
— 1/2, the distributions display a marked asymmetry. Also apparent, is a pronounced shifting of 
the high density peak to successively lower densities as A/77 is increased. By contrast, however, 
the position of the low density peak is relatively unaffected by changes in A/77. Measurements 
of the critical density (calculable simply as the first moment of these distributions), agree well 
with the exact values obtained from equation §. 

That the asymmetry of the distributions vanishes for A/77 = — 1/2, is demonstrated in 
figure |] where a comparison is shown between the critical density distribution and the mag- 
netisation distribution of the 2-D spin-- Ising model at its exact critical point, as previously 
obtained by Nicolaides and Bruce Jl7|- The density data has been expressed in terms of the 
scaling variable x = aM~ 1 L^' v [p — p c ], and the data collapse has been effected by a single 
scaling (choice of the non-universal scale factor clm) such that both distributions have unit 
variance. Clearly, the overall quality of the data collapse is impressive, except for small depar- 
tures near the vestiges of the density distribution. These however, can be traced to the fact 
that the L = 20 system is not quite sufficiently large to allow proper sampling of the high and 
low density tails of the density distribution. For this reason, detailed analysis of the corrections 
to the density distributions was carried out for the L = 30 system where this problem is much 
less evident. 

As previously noted, the antisymmetric component of the density distributions may be 
isolated simply by antisymmeterising the density distribution about p c , whose value is obtain- 
able from equation |[ The results of applying this procedure are depicted in figure |5] for the 
L = 30 system at A/77 = —5/12 and A/77 = 0- m the former case, the size of the antisymmetric 
component is quite small, while in the latter it is considerably larger. These antisymmetric 
components are replotted as functions of the scaling variable x in figure ^|, together with the 
predicted universal form of the field mixing correction, prescribed in equation [15] and obtained 
from Ising model studies [|17], ||. The data have all been brought into coincidence by single 
scalings of the non-universal scale factor aju and the ordinate. In both instances the antisym- 
metric component of the density distribution maps extremely well onto the predicted universal 
form. 

Turning now to measurements of the significant mixing parameter, the value of s can be 
estimated according to the procedure outlined in the previous section. We shall describe this 
procedure for the case A/77 = —5/12, for which the antisymmetric component of the density 
distribution is rather small. Some additional complications arise when the antisymmetric 
component is larger, and these are discussed separately in the context of the A/77 = data. 

In order to measure s , both the energy function < u(p) > —u c and the density distribution 
Pl(p) are required. The energy function is displayed in figures [j] for the case A/77 = —5/12. Also 
included in figure [7|is the function r(p — p c ), where r (cf equation |) is the limiting critical slope 
of the coexistence curve in the plane of 77 and £, whose value (for a particular choice of A/77) 
is calculable from equation [|| For the present case, A/77 = —5/12, one finds r = —1.054092. 
Now, according to equation [14], s is simply the single scale factor required to match the 
function Ap L (p) = —s-^- {pl(p) [< u(p) > —u c — r(p — p c )]} to the antisymmetric component 
of the density distribution. To the extent that O(s) correction to the energy function and the 
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density distribution can be neglected, the function Apx(p) will itself be antisymmetric and 
thus comparison of the scales of the two functions can be effected directly, notwithstanding 
the relative enhancement of statistical uncertainties in Api(p) that results from taking the 
numerical derivative. Figure [8] shows both Api(p) and the measured antisymmetric component 
of the density distribution. The two data sets have been brought into correspondence by a 
single scaling (choice of the ordinate) implying a value s = —0.23(2). This estimate is in quite 
good agreement with the exact value of —0.204.... 

Finally, we address briefly the situation where the antisymmetric component of the density 
distribution is not small. Under these circumstances, O(s) terms make a significant contri- 
bution to both Pl(p) and u(p). These terms (manifesting themselves as 0(s 2 ) corrections to 
Api(p)) prevent a direct comparison of Ap^(p) with the measured antisymmetric component 
since Api{p) will no longer itself be antisymmetric in p — p c . Notwithstanding this, an esti- 
mate of s can be obtained by the simple expedient of antisymmeterising Ap L (p) about p c , and 
comparing the scales of the two functions. Carrying out this procedure for the X/t] = data 
yields an estimate s = —0.38(1) to be compared with the exact value s = —0.331.... The bulk 
of this discrepancy, some 15%, is presumably attributable to the neglect of field mixing terms 
of order s 2 and higher, which are not included in equation 115. 



4 Conclusions 

The simulation results presented in this paper confirm that in systems with broken particle- 
hole symmetry, the field mixing phenomenon gives rise to an antisymmetric correction to the 
limiting form of the critical density distribution. An accurate determination of the functional 
form of this correction has demonstrated it to be in excellent quantitative agreement with the 
prediction of the mixed- field finite-size-scaling theory of references [7], |8|]. As such the results 
constitute substantial corroboration of that theory. 

A method for estimating the significant mixing parameter s was also assessed, and found 
to yield accurate results, provided the overall magnitude of the field mixing correction is not 
too large. In view of the finding || that the size of the field mixing correction in the Lennard- 
Jones fluid is indeed rather small (constituting only approximately 5% by weight of the density 
distribution for a system of 400 particles), the method described should in principle permit 
quite accurate evaluation of s in realistic (off-lattice) fluid models. Of course (and as noted 
in the introduction), the feasibility of such studies is contingent upon an ability to first locate 
accurately the critical point for the model of interest. However, as was demonstrated in refer- 
ence ||, this can also be achieved via a finite-size-scaling analysis of the density distribution, 
provided one works within an ensemble (such as the grand canonical ensemble) that affords 
adequate sampling of density fluctuations. 

With these points in mind, it would be of considerable interest to extend the present 
programme of simulation studies to an investigation of the microscopic factors influencing the 
degree of field mixing in realistic fluid models. Recent theories have predicted that the size 
of the significant mixing parameter is strongly influenced by three-body forces ||, 0. Such 
interactions could certainly (albeit at greater computational expense) be incorporated into 
off-lattice fluid simulations, allowing their role in the field mixing process to be assessed. 
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Figure 1: The critical density p c of the asymmetric lattice gas model expressed as a function 
of the coupling ratio A/77 



Figure 2: The significant field mixing parameter 's' expressed as a function of the coupling 
ratio A/77 



Figure 3: The L = 20 critical density distribution for various values of the coupling ratio A/77. 
Statistical uncertainties are smaller than symbol sizes; the curves are simply guides to the eye. 
All distributions are normalised to unit integrated weight. 
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Figure 4: The critical density distribution for the L = 20 system at A/77 = — 1/2, expressed 
as a function of the scaling variable x = L /3 / 1/ a^/(p — p c ). Also shown (full curve) is the 
critical magnetisation distribution of the 2-D Ising model obtained in ||17|| . The non- universal 
scale factor implicit in the definition of the scaling variable, has been chosen so that both 
distributions have unit variance. 



Figure 5: The measured antisymmetric component of the L = 30 critical density distributions 
at A/77 = —5/12 and A/77 = 0. The lines are guides to the eye. 



Figure 6: The data of figure |5] re-expressed in terms of the scaling variable x = LPI u a N }(p — p c ) 
and shown as the data points. The solid curve represents the prediction following from equation 



15, utilising predetermined Ising forms 17 ]. The measured correction data have been brought 



into coincidence with the predicted form via appropriate choices of the non-universal scale 
factor cim and the ordinate. 



Figure 7: The measured form of the critical energy function < u(p ) > — u c for the L = 30 
system at A/77 = —5/12. Statistical uncertainties do not exceed the symbol sizes. Also shown 
(solid line) is the function r(p — p c ), with r = —1.054092. 



Figure 8: The measured function Ap L (p) = —s-§^ {pl(p) [< u(p) > —u c — r(p — p c )}} for the 
L = 30 system at A/77 = —5/12 and shown as crosses (x). Also shown as circles (o) is 
the measured antisymmetric component of the critical density distribution, cf. figure |^. The 
matching shown was effected by a choice of the ordinate implying s = —0.23(2). 
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